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摘 要 : 为 了 得 到 一 类 中 立 型 灰色 随机 分 布 时 灌 系 统 的 指数 鲁 棒 稳 定性 ， 本 文 利用 Lyapunov- 
Krasovskii 泛 函 法 、 灰 矩阵 的 连续 矩阵 覆盖 的 分 解 技 术 和 HI6 人 公式， 分 别 得 到 了 以 非 线性 和 矩 
阵 不 等 式 和 线性 矩阵 不 等 式 (LMI) 表示 的 该 系统 指数 鲁 棒 稳 定 的 时 澡 依 赖 性 判 据 。 对 非 线性 矩 
阵 不 等 式 判 据 ， 我 们 给 出 了 一 般 性 算法 ， 解 决 了 非 线性 矩阵 不 等 式 判 据 不 便于 实际 应 用 的 问题 。 
数值 例子 表明 ， 本 文 所 给 判 据 是 有 效 的 ， 且 系统 的 指数 稳定 性 和 时 滞 ， 随 着 绝对 灰 度 矩阵 的 谱 范 


数 的 增 大 而 减 小 。 
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1 引言 


由 于 随机 微分 系统 在 自然 、 社 会 和 科技 等 领域 具有 广泛 的 应 用 ， 所 以 ， 近 三 十 多 年 
来 ， 关 于 随机 微分 系统 的 稳定 与 控制 问题 ， 一 直 是 很 多 学 者 关注 的 焦点 问题 ， 并 取得 了 
许多 有 价值 的 成 果 上 看。 其 间 ， 在 二 十 世纪 八 十 年 代 ，Kolmanovaskii 和 Nosov 基于 化 学 工程 
和 航空 理论 发 展 的 需要 ， 还 建立 了 一 类 中 立 型 随机 泛 函 微分 系统 ， 并 研究 了 该 系统 解 的 
存在 性 、 唯 一 性 饭 、 稳 定性 和 渐 近 稳定 性 四 的 问题 。 此 后 ， 一 些 学 者 又 利用 Lyaponuv 泛 函 
和 Razumikhin 技术 ， 研 究 了 中 立 型 随机 系统 的 指数 稳定 问题 ， 得 到 了 一 些 指数 稳定 的 代数 判 
lte], 20074F, Randjelovic 和 Jankovic 则 基于 一 个 关于 向 量 范 数 的 产 阶 矩 不 等 式 ， 给 出 了 中 
立 型 随机 系统 的 广 阶 矩 指数 稳定 的 代数 判 据 四 ， 且 该 判 据 还 是 文献 [8] 的 结果 的 推广 。 然 而 ， 
关于 中 立 型 不 确定 随机 系统 的 研究 ， 目 前 仅 有 少量 的 报道 。 如 文献 9] 在 建立 了 中 立 型 随机 
系统 的 LaSalle 不 变 原理 的 基础 上 ， 应 用 该 原理 讨论 了 中 立 型 不 确定 随机 系统 的 随机 渐 近 稳 
定性 和 几乎 必然 指数 稳定 的 代数 判 据 。 文 献 [10] 则 给 出 了 中 立 型 不 确定 随机 系统 的 鲁 棒 镇 定 
和 五 控制 的 一 种 状态 反馈 控制 器 的 设计 方法 。 至 于 有 关中 立 型 随机 分 布 时 滞 系 统 的 研究 ， 晶 
前 还 没有 相关 报道 。 灰 色 系 统 理 论 是 从 另 一 个 角度 来 研究 和 处 理 自然 和 社会 系统 中 的 不 确定 问 
题 的 一 个 典型 工具 ， 它 能 有 效 处 理 不 确定 系统 的 稳定 与 控制 问题 。 但 是 ， 中 立 型 灰色 随机 分 布 
时 滞 系 统 的 鲁 棒 稳 定 与 控制 问题 ， 还 没有 人 涉及 。 另 外 ， 在 研究 中 立 型 不 确定 非 随 机 系统 时 ， 
所 采用 的 一 些 模型 变换 方法 ， 不 适用 于 处 理 中 立 型 不 确定 随机 系统 的 问题 ， 因 为 中 立 型 随机 系 
统 的 解 ， 一 般 是 不 可 微 。 因 此 ， 本 文 将 用 Lyapunov-Krasovskii 泛 函 法 、 灰 矩阵 的 连续 矩阵 履 
盖 的 分 解 技术 和 一 些 基 本 不 等 式 ， 研 究 一 类 中 立 型 灰色 随机 分 布 时 滞 系 统 的 指数 鲁 棒 稳 定性 问 
题 ， 给 出 其 指数 鲁 棒 稳 定性 的 线性 矩阵 不 等 式 条 件 。 
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2 系统 描述 


设 T>0,C(-7,0];R") 是 定义 在 [-7,0] 上 并 取 值 于 R* 上 的 连续 函数 o 的 全 体 。(Q, Z, 
{F}i>0,P) 是 一 个 具有 自然 流 {天 }i>o 的 完备 概率 空间 ，w(t) 是 定义 在 概率 空间 上 的 一 维 的 
标准 布朗 运动 。L.([ 一 7,0]; Rz) 表示 天 可 测 的 取 值 于 C([--7,0];R"*) 上 的 随机 变量 & = {€(8) : 
-r < 0 < 0} 的 全 体 且 使 得 


sup E|€(8)|? < co. 
—r<0<0 


zt = {x(t + 0) : -r < 0 < 0} ERT C([-7, 0]; R”) 上 的 随机 过 程 。 令 | :| 和 外. 半分 别 表示 一 
个 向 量 和 一 个 矩阵 的 欧 氏 范 数 。 如 果 4 是 一 个 对 称 矩 阵 ， 用 符号 4< 0 和 4 < 0 分 别 表示 4 是 
一 个 半 负 定 和 负 定 矩阵 ， 并 用 Amax(4) 和 Amin(A) 分 别 表示 它 的 最 大 和 最 小 特征 值 。 另 外 ， 
& AT 表示 4 的 转 置 ，@，4(@), SAAG) 分 别 表示 一 个 灰 数 、 灰 和 矩阵、 白化 数 和 白化 矩阵 。 
定义 1 如 果 灰 矩阵 4(@) 中 的 灰 元 素 的 数值 覆盖 全 是 连续 覆盖 ， 那 么 这 个 灰 和 矩阵 A(@) 的 
矩阵 覆盖 被 称 为 连续 矩阵 覆盖 Q0， 并 称 这 样 的 灰 矩 阵 为 具有 连续 覆盖 的 灰 矩 阵 或 简称 区 间 灰 


矩阵 。 
基于 下 面 的 一 个 中 立 型 随机 分 布 时 滞 系 统 


d|z(t) — Fz(t — r)] = [Az(t) + Bz(t — r) + f, G(t — s)z(s)ds]dt 
+[Cz(t) + De(t —7)|dw(t), t> 0, (1) 
z(t) =E(t), €€ LE ([-z, OJ R*), —+r <t<0. 


我 们 给 出 中 立 型 灰色 随机 分 布 时 湾 系 统 及 其 指数 鲁 棒 稳定 性 的 定义 。 

定义 2 如果 系统 (1) HRSG F, A, B, C, D 和 权重 矩阵 G(t) 中 至 少 有 一 个 灰 和 矩阵 ， 那 
么 系统 (1) 被 称 之 为 中 立 型 灰色 随机 分 布 时 滞 系 统 。 

本 文 考虑 下 述 一 个 中 立 型 灰色 随机 分 布 时 灌 系 统 


d[z(t) ~ F(®)x(t — r)] = [A(®)z(t) + B(@)a(t — r) + e G(t ~ s)x(s)ds|dt 
+[C(®)zx(t) + D(@)x(t — 7)]dw(t), t>0, (2) 
a(t)=€(t), Ee L%,([-7,0];R"), -7r<t<0, 


其 中 z(t) ER", A(®), B(@), C(®), D(@) 和 F(@) E: n x n H: LAE EAE BE 3: AI, 
它们 的 连续 矩阵 覆盖 或 白化 矩阵 集 分 别 为 


Ar = [La, Ua] = {A(®) = (Bais)nxn : G; < Sais < Gj, 4.7 =1,--- ,n}, 
Bi = [Ls, U,] = {B(Q) = (viz) nxn : bij < Ooij S bijs 4,7 =1,--- ny, 
Cy = [Le, Ue] = {C(®) = (@eij nxn : cij < ois S Gj, bj = 1 ,n}, 
Dr = (La, Ua] = {D(®) = (@aij nxn :di; < Baij < diz, if = 1,::: n}, 
Fy = [L/,U/] = {F (8) = (Šfij)nxn fy < yey < Fiz, Lj=1, np. 


G(t) 是 一 个 nxn A RL. ENRE x(t) = E(t) e L} ([—7,0];R") F, ABE (2) HY 
fF a(t; €) 相应 于 初始 值 E(t) = 0， 其 平凡 解 z(t;0) = 0. 
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定义 3 ”对 于 任意 的 白化 矩阵 4(@) € Ar, B(S) € Br, C(®) € Cr, D(&) € Dr AM F(&) € 
Fr 及 任意 的 & e LZ ([-7,0]; R")， 如 果 存 在 正常 数 ro 和 天 ， 使 得 系统 (2) 的 解 满足 


Elz(t;é)| < Ke™ sup Elé(O)|, t20, 
—T<0<0 
或 等 价 地 
lim sup [SinE |z (t; £) P] < —ro, 
t—oo t i 
那么 称 系统 (2) 是 均 方 指数 鲁 棒 稳 定 的 ， 如 果 有 


lim sup | 5 In |z(t; 8l] <——, as., 
t—oo 


那么 称 系统 (2) 是 几乎 必然 指数 鲁 棒 稳 定 的 。 
512197] 如 果 W(@) 是 一 个 m xz 的 灰 矩 阵 ， 它 具有 连续 矩阵 覆盖 


Wy = [Lw, Uw] = {W(®) = (@wiz)mxn : Wij S Quiz < Wiz, 1=1,---,m, = 
那么 对 任意 的 白化 矩阵 W(@) e Wr 都 有 
W(®)=Ls+AW HB. |AWI < |JU, - Lull 
其 中 
Uw = (Wij)mxn, Lw = (Wij)mxn AW = (ky)mxn, kij = Wig — Wiz, Vag € (0,1), 


(0, 1] FEAL ARK y(@) 的 连续 数值 覆盖 。 
引 理 2 (Schur 补 ) 对 于 给 定 的 对 称 和 矩阵 


g= Sir S12 | 
S21 S22 
其 中 SeR"X", S, ER”, r<n, 有 


S<06Sy<0, Sn- Ssa Ssh < 0. 


3 ”主要 结果 


定理 1 设 k 二 上 Lf 川 + Us 一 万 | < 1。 如 果 存 在 正常 数 e1,… ,eo REM Q fü R, HH 
得 


Amin(R) > (ez + Eg Ds p IGO, (3) 


Mex Mı + piln La + LT La = LIL; <0, 
LT + LẸ Le = LF La M2 + Haln 
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HH In ER n x nn 单位 矩阵 

Mı = La + LT + e3L7 La + (1+e7)L7L.+Q4+7R, 

Hy = £1 + €2 + 2||U, — Lall + (gz! + €5) [Ua — Lall? 
+ 2||Ue > Lelll|Lcll + (1 + €g* + €9) Ue — Lell?, 

Mz = (€4 + €6)L} Ly + (1 +eg)LG La — LF Ly — LẸ Lf — Q, 

H2 = E7 (Us ~ Lell? + (gs! + e5 Uy — Lill? + €6 (Oy — Ly ll? +210; — Lylli| L|) 
+(1+e7! +es Ua — Lall? + 2|JU; — LyllliLol + 2U» — LollllLsll 
+ 205 ~ LyllllU, — Loll + 2a — LallllLall, 


那么 系统 (2) 是 均 方 指 数 鲁 棒 稳定 的 ， 且 Lyapunov 指数 为 -ro。 而 ro 满足 ; ro > 0, ket < 
1 及 
ro(1 + k?) + Amax(M) + roAmax(Q)re"” + roAmax(R)r2e"” < 0. 


证 明 对 任意 固定 的 E e LZ ([-7,0];R") 和 白化 矩阵 4(@) e Ar, B(&) € By, C(B) € 
Cr, D(®) € D; K F(@) € Fr, WHE (2) HR r(t; £) = z(t), Hid 


t 
y(t) = z(t —7), z(t) = / G(t — s)x(s)ds. 
t-r 
MF (t,£) € R, x R”, W Lyapunov-Karsovskii ZZ % 
V(t, z(t)) = [z(t) - FEW] [z(t) — F@)v)] 
t 0 pt 
x" (s)\Qz(s T 
+ f (s)\Qz(s)ds + f. I: 2Z (s)Rz(s)dsd0, 
则 利用 It6 公式 中 的 弱 无 穷 小 算 子 LV， 沿 着 系统 (2) 的 白化 系统 的 迹 计 算得 
LV (t,2(0) = TOQ ETRE) -TOU — #78) Ra(s)ds +22" (0) A(@)z(9) 
—2y7 (t)F7 (8)A(B)TE) + 2x7 (] B(@)u(t) + 2x7 (t)z(t) 
—2y7 (t)F7(&)B(®)y(t) — 2yT (0 F(@)z(t) + zT (t)CT(&)C(&)=(t) 
+2y? (t) DT (@)C(®)x(t) + y” (t) DT (®)D(®)y(t). (4) 
于 是 ， 由 引 理 1 可 推 得 
227 (t)A(®)a(t) = x7 (t)(La + LT) x(t) + 227 (t A Az(t) 
< z” (t)(La + La)z(t) + 2||Ua — L, zT (t)x(t), (5) 


227 (t)B(®)y(t) < 27 (t)Ley(t) + y7 (tÜ) Ë z(t) + e127 (t)z (D) 


+e; Us — Lell yT (Du(t), (6) 
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2z7 (t)z(t) < e2xT (t)x(t) + ex zT (t)z(t), (7) 


~2y™ (t) FT (@)A(@)z(t) = —y (t)LF Lax(t) — x7 (LT Lsy(t) — 2y7 (t) A FT Laz(t) 


— 2y7 (t)L? AAx(t) — 2y7 (AF? AAz(), (8) 
~2y" (t) AF? Laz(t) < e327 (ü) LT Laz(t) + es Us — Lfl uT yl), (9) 
-2y (Lf AAx(t) < e4l|Ua — Lall?x" (t)x(t) + eg ty” LF L ry(t), (10) 
—2y? (t)AF? A Az(t) < €5||Ua — Lall?x7 (t)z(t) + es ` lU; — Lely (y(t). (11) 


类 似 地 ， 
—2y7(t)F7(@)B(®)y(t) < —y (t)(LF Ly + L5 Ls)y(t) + (2U; ~ LsllllLel 
+2||U, — LalillL sll + Us — Lsls — Lolly y(t), 2) 

—2y7 (t)F(@)z(t) < eey” (t) T (@)F(@) y(t) + =s '27 (t)z(t), (13) 
YTF (SFO) Sy (LF Ly y(t) + (2lJU; — Lele + Ws — Lell?)y7 @y@), (14) 
zT (CT (&)C(ë)=(t) < z! (t) LZ Lez(t) + (2||Ue — Lell||Lel] + Ue — Lell?) (ald), (15) 
2y (t) DT (&)C(ë)z(t) < yY (t) LT Lex(t) + x? (E) LT Lay(t) + erzT (t) LT Lez(t) 

+7 Ua — Lally” (yd) + esy” (t)La Laylt) 

+es ||Ue — Leh 2T (t)x(t) + eollUe ~ Lell?27 (t)z (t) 

+ €9*|\Ua — Lall’ y” (t)u(t), (16) 


Y” (HDT (&)D(@)y(t) < y” (La Lay(t) + (2\|Ua — LallllLal + Ua — Lal?) (Hy(t). (17) 


另外 ， 
_ f ü zT(s)Rz(s)ds < —Arin(R) Ë zT (a)=(s)ds, as 
及 由 Halder 不 等 式 知 
27 (t)a(t) < r J i IG(t~ s)z(s)?ds < r sup 1G(s)P S Iz(s)|2ds, T 


所 以 ， 将 式 (5)-(19) RAR (4) 中 整理 后 ， 再 注意 到 M, Mi, Mo, p, woe 的 定义 及 条 件 (3) 可 得 


x(t) 


LV (t,x(t)) < (x7 (t), yT (t))M 
y(t) 


| < Amax(M)[|x(t)|? + |u(t)|2]. (20) 
FES w(t, s(t) = erotyY， 则 由 It6 公式 知 


Ew(t,2(t)) = Ev(0,2(0)) + J Bro (s, zed (21) 
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其 中 
Ly l(t, z(t)) = roe" V (t, x(t)) +e” LV (t, x(t). 


而 根据 Cauchy 不 等 式 和 引 理 1 知 
|z() — F(®)| < (1 + k2) [la(e) |? + y(t)? J, 


且 注意 到 ro(1 + k?) + Amax(M) < 0， 所 以 易 得 


L(t, 2(t)) < e { [ro(1 +k) + Amax(M)]lz(bP+roAmax(G) f Iz(s)Pds 


0 t 
Amax(R 2dsd0 $. 22 
+rodman(R) Í Í Ie(9easao) (22) 
# C = sup_,<e<o ElE(0)|2?， 则 可 推 得 
Ey (0, x(0)) < [2(1 + k?) + TAmax(Q) + T?Amax(R) JC, (23) 
[om f° Bleavas < rew | eBla(s)as + eo, (24) 
0 


f = Ji I E|z(u)|?dud6ds < + fe 7 Bla (u)|2duds 


< ° f erosE|z(s)|?2ds + r3ero7C. (25) 
0 
于 是 ， 由 式 (22)、(24) 和 (25) 可 得 


t 
f ELy(s,2(s))ds < [ro(1 + k?) + Amax(M) + roTÀmax(Q)e"%” + 7977 Amax(R)e”?” | 
0 


x | erosElz(s)| ?ds + roXmax(Q)r2erorC + roAmax(R)T3e"7 C. (26) 
0 
此 外 ， 利 用 基本 不 等 式 2ab < ca? + g 102 (e > 0)， 并 注意 到 天 的 定义 ， 我 们 有 
Ey(t, z(t)) > Efe lelt) — F(@)u(t)| ] 
> eret[(1 — k)E|x(t)|? — (k — DE|F(ë)u(#)| ] 
er [ (1 — k)E|z(t)|? ~ (1 — k)KE|y(t)|? ]. (27) 


IV 


这 样 ， 结 合式 (21). (23). (26) 和 (27) É ro 的 定义 便 有 
(1 — k)e”*E|z(t)|? < Co + (1 — k)ker°tE|y(t)|2 


< Co+(1—k)ke"* sup Elz(t +0), (28) 
—+r<0<0 


其 中 
3 707 


Co = [2 (1 T k?) T TÀmax(Q) + T2Àmax(R) T ToAmax(Q)T7e ror p roAmax (R)T e ]C. 
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对 于 任意 的 0 < t< T, T> 0， 由 式 (28) 知 


erotg|z(t)|2 < Co +k sup [ert sup E|z(t+ o)? | 
l-k O<t<T —7<6<0 
= ow +k sup | sup erm(s-9)E |z(s)| | 
l-k o<t<T | t-r<s<t 
< = +ke7 sup [e’E |x(t)|? ] s=t+6. (29) 
1—k —r<t<T 
显然 ， 对 于 一 + <t <0, (29) 也 成 立 。 于 是 
sup [e’*E|z(t)|?] < Co + ke"? sup [e”*E|z(t)|?]. (30) 
—r<t<T l-k —r<t<T 


这 样 ， 由 式 (30) ro 满足 的 条 件 可 知 


(31) 


Ele(t)? < 2 


e rot 
(1 — k)(1 — kero7) š 
即 系统 (2) 是 均 方 指数 鲁 棒 稳 定 的 。 

注 1 在 实际 应 用 时 ， 要 确定 定理 1 中 的 判断 矩阵 M 是 否 负 定 ， 关 键 是 能 否 找到 使 它 负 定 
的 参数 s; (i =1,2,--- ,9) 和 正定 矩阵 Q@ 及 有 R， 而 用 手工 调试 寻找 的 办 法 显然 是 不 易 的 。 不 过 ， 
我 们 可 以 先 让 

Q=kiln, ki>0, R=keln, ke > (ez'+e5')7 sup ||G(s)|l’, 
0<s<r 


这 样 ， 便 可 以 把 此 问题 看 作 是 一 个 关于 e;，k1 和 ko 的 多 元 非 线 性 矩阵 函数 的 负 定 问 题 。 昌 然 此 
时 仍 不 易 有 效 地 采用 手工 方式 获取 一 组 e;,，k1 和 kz 值 ， 但 本 文 可 以 采用 类 似 于 文献 [12] 中 的 自 
由 参数 的 优选 方法 ， 通 过 编制 一 个 程序 ， 让 计算 机 自动 搜索 一 组 使 判断 矩阵 负 定 且 具有 较 快 收 
SER en ki 和 ks 值 。 当 然 ， 也 可 以 先 取 定 有 和 ko2， 然 后 编程 选取 自由 参数 e;。 自 由 参数 的 
优选 算法 为 : 

步骤 1 输入 已 知 矩 阵 和 常数 ; 

步骤 2 给 定 最 大 搜索 计算 次 数 ， 如 i = m。 选 用 一 个 for 循环 指令 ; 

步骤 3 ”让 计算 机 随机 产生 一 个 1 维 向 量 e (1 表示 非 线 性 判别 矩阵 中 的 自由 参数 的 个 数 ，e 中 
的 1 个 随机 数 是 区 间 [0,1] 上 的 数 ); 

步骤 4 计算 判别 矩阵 M 及 其 最 大 特征 值 入 ; 

步骤 5 如 果 和 和 <0， 让 入 = 入, ct =e, A, BAAR: 

步骤 6 如 果 i 二 1, 令 max 人 入 * = 入，maxe* = e*， Ail, $ maxà** = A*, maxe** 一 Er; 

步骤 7 如 果 max 入 * < max 入 ， 令 max 入 = max à**, maxe* = maxe™*, BM, KEY 
UR 2; 

步骤 8 Hi=m, AREA., 否则 ， 返 回 步骤 2。 

下 面 ， 我 们 基于 定理 1， 推 出 一 个 由 线性 矩阵 不 等 式 (LMI) 表示 的 指数 鲁 棒 稳 定性 判 据 。 

推论 1 WR = ||Ly|| 十 上 VF 一 Lil| < 1。 如 果 存 在 正常 数 al,… , a9 和 正定 矩阵 Q@ 使 得 矩阵 


Ni La + LT La 一 LTL; Q 
N := LT + LT L. 一 LẸ La N> + Ns 0 < 0, 
QT 0 —J 
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其 中 
Ni = La + L7 + Li Le + Q 站 (2||Ua š; Lall T 2||U. tk LecllllLell + lU; = Dal? Mn, 
Na = Lå La — LẸ Ly — Lp Ly — Q + [||Ua — Lall? + 2U — Lyl\l|Loll + 2U — LollllL sll 
+2|JU; — Lillis — Loll + 2||Ua — LallllLall Zn, 
Ns = [oils — Loll? + (a3 + os)||Uy — Ly ll? + a6(|Up — Esl? + 2U; — LENEA) 
+ (a7 + ag)||[Ua 一 Lall? Zn + (a4 + ag) LF Ly + as LI La, 
而 
J =diag(aiIn, a2In, azn, osl,, oaln, O5In, oeln, 7In, A8In, ooln) 
FER kE E, 


Q = (In ool, T Í sup |G(s)|2I,, LZ, Ua — LallIn; 
OLST 
Ua — Lallin, 7, Í sup ||G(ə)|T,, LZ, [Ue — LellIn, lUa ~ Lallin) 
O<s<r 


是 一 个 单行 块 矩阵 ， 那 么 系统 (2) 是 均 方 指数 鲁 棒 稳 定 的 ， 且 Lyapunov 指数 为 -ro。 而 ro 满 
Æ: ro > 0, kT < 1 及 


ro(1 +k?) + Amax(N) + roAmax(Q) re" + ror%e"" (az! +ag') sup ||G(s)||? < 0. 
O<s<r 


证 明 BOTH 1 PH 
R= (eg! +e51)7 sup |G(s)|21,, a:=ez', 1=1,3,5,7,9, oj =e), J =2,4,6,8, 
O<s<r 
同时 ， 根 据 引 理 2 进行 变换 ， 即 可 由 定理 1 证 得 推论 1。 
定理 2 ”如 果 定 理 1 或 推论 1 中 的 条 件 成 立 ， 那 么 系统 (2) 又 是 几乎 必然 指数 鲁 棒 稳 定 的 。 


证 明 若 定 理 1 或 推论 1 中 的 条 件 成 立 ， 即 系统 (2) 是 均 方 指数 鲁 棒 稳定 的 ， 则 存在 正常 
Bro Al K 114 


Bla(t;6)? S K sup Eé(0)2e”*  €ELF,([-7,0];R"), t>0. 
—r<0<0 
现 记 
b= (Lp + Uy -i <1 #= min {ro, T7}, In I 21, 


则 对 于 任意 的 5 © (0,7) 和 任意 的 整数 k > 1， 应 用 著名 的 Doob WERTERA Cauchy PEA 
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P(w : sup jz (kr +s)— F(@)zkr+. |° 5 er) 
O<s<r 


E|x(kr +7) 一 F(@)zk+-++|° < ü + ||F(®)|l?) (Ela(kr + 7)|? + E|zk+ ++ |?) 
< - 
= e- (F—ó)kr e—(f—ó)k+ 


(1 + 2) [He-reter+7) sup E|€(6)|? + Ke ro(kr+r+0) sup EJE(6)1?] 
—r<6<0 一 r<6<0 
e—(f—ó)kr 


= K(1 41?) (e777 + e-ro(r+0))e(f—ro)kr—ëkr sup E\€(0)|? 
一 r<9<0 
< 玫 (1+D(L+ernr)e- or sup El|é£(6)Ë, (32) 
一 T<O<0 


进而 ， 由 著名 的 Borel-Cantelli 引 理 知 ， 对 于 几乎 所 有 的 we Q 和 除了 有 限 个 正 整数 大 之 外 的 所 
有 正 整 数 ， 都 有 
sup |a(kr +s) — F(®)zer+s|° < ek (33) 
O<s<r 
因此 ， 存 在 一 个 正 整数 fo(w)， 使 得 对 于 几乎 所 有 的 we Q， 包 括 尸 零 可 测 集 ， 当 玉 > kh, 
A (33) 都 成 立 。 也 就 是 说 ， 如 果 上 > por， 那 么 ， 对 于 几乎 所 有 的 we 2， 都 有 
jz( — FP(ë)=,|° < e F-8)E, 


HF |z (t) —F(@)z,|2 在 区 间 [0, kor) 上 是 有 限 的 ， 所 以 ， 存 在 一 个 有 限 的 正 数 C* = C (w), 
使 得 对 于 几乎 所 有 的 we Q t> 0， 都 有 


|x(t) — F(@)z,|? < Cte“ 0)t. (34) 


男 外 ， 基 于 式 (34)， 类 似 于 式 (30) 的 推导 过 程 ， 我 们 可 以 得 到 : 对 于 任何 的 实数 全 > 0， 都 有 


(1 ~ le*-5)7) Pe [e-t a(t)? ] < =. 
进而 ， 得 a 
(F—)t 2 
°" POMS Gaya ey)’ 
这 样 便 有 


, 1 l (f — ó) 
Hm sup [zm |z (8; 6)|| < —— Gs 
Fë, #4550, WA 
lim sup Ë In |z(t; 8l] < -5 a.8., 
too 
这 就 说 明了 定理 2 是 成 立 的 。 


注 2 事实 上 ， 在 上 述 定理 或 推论 的 判 据 中 ，||Ui — Lill, i= a,b,c,d, f TRACERS (2) 的 
灰 参 数 的 绝对 灰 度 构成 的 矩阵 的 范 数 。 如 果 这 些 灰 参 数 的 绝对 灰 度 是 变化 的 ， 那 么 ， 我 们 可 以 
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基于 推论 1 的 判 据 ， 来 探讨 灰 参 数 的 绝对 灰 度 的 变化 对 系统 的 指数 稳定 性 和 时 滞 长 度 的 影响 情 
况 。 文 中 把 灰色 和 矩阵 的 连续 覆盖 进行 如 此 分 解 的 好 处 就 是 ， 既 可 以 借 此 讨论 绝对 灰 度 矩阵 的 变 
化 对 系统 的 指数 稳定 性 和 时 滞 长 度 的 影响 情况 ， 又 可 以 减少 或 避免 灰 矩 阵 中 所 包含 的 信息 量 的 
损失 或 破坏 。 

注 3 如 果 令 系统 (2) 中 的 G(s) = 05k F(@) = 0， 我 们 可 以 利用 定理 1 或 推论 1， 分 别 得 到 
与 系统 (2) 相对 应 的 中 立 型 灰色 随机 时 滞 系 统 和 灰色 随机 分 布 时 滞 系 统 的 指数 鲁 棒 稳 定 的 判 
据 ， 这 里 不 再 详 述 。 而 对 这 两 类 灰色 系统 的 指数 鲁 棒 稳 定性 的 研究 ， 目 前 也 没有 相关 报道 。 


4 ”数值 例子 


例 1 考虑 一 个 二 维 的 中 立 型 灰色 随机 分 布 时 光线 性 系统 ， 其 灰 参 数 矩 阵 的 连续 矩阵 覆盖 
Æ Fi, Ar, Br, Cr ® D: WE FREE RABE RAD AA 


—0.285 0.113 —0.275 0.125 一 3.870 —0.531 

Ls = ; Us > ; La = 5 
0.112 —0.283 0.123 —0.273 0.513 —3.760 
—0.639 0.138 —0.629 0.146 

L, = 》 U, = 5 
0.106 —0.636 0.117 —0.626 


0.136 —0.113 0.146 —0.104 0.417 —0.169 
Le = , Ue = , La = ’ 
—0.103 0.125 —0.092 0.136 —0.112 0.422 


系统 的 时 滞 长 度 为 = 0.1. 
显然 ， su G(s) =e. FOZHI NABH 


Q=kiln, Amin(R) = (eg*+e5")7 sup |G(s)|, 
O<s<r 

那么 ， 利 用 前 面 所 给 的 算法 ， 通 过 对 M < 04 Matlab Editor/Debugger 环境 中 编程 计算 ， 获 
得 的 使 判别 矩阵 M 较 快 收敛 的 参数 esi， i= 1,2,… 9, k 及 由 此 而 得 的 最 大 的 负 特 征 值 和 指 
数 收 敛 速率 分 别 为 

gi = 0.0779, £2 = 0.199, £3 = 0.0076, sd4 = 0.4887， &5 = 0.0225, 

g é = 0.1090, £7 =0.1747, sg = 0.0742, <£ = 0.7983, k, = 3.0000, 

Amas (M) = —1.2372, ro = 0.7685. 

从 上 述 的 数值 结果 可 以 看 到 ， 根 据 定理 1， 我 们 可 以 判断 出 本 例 所 描述 的 系统 是 指数 鲁 棒 稳 


定 的 ， 且 有 较 快 的 指数 收敛 速率 。 也 就 是 说 ， 定 理 1 所 提供 的 判 据 都 是 有 效 的。 当然 也 说 明了 


推论 1 的 判 据 的 有 效 性 。 
下 面 ， 我 们 基于 推论 1， 利 用 例 2 来 说 明 灰 参数 矩阵 的 绝对 灰 度 矩阵 的 变化 对 系统 的 指数 稳 


EPERRA KEKE , 
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例 2 ”假设 在 例 1 的 二 维系 统 中 ， 系 统 的 灰 参 数 矩 阵 的 绝对 灰 度 矩阵 都 是 变化 的 ， 且 绝对 灰 
度 矩 阵 满足 lz — Lill <6, ó > 0, L, #1G(s) 同 例 1, ¿= a,b,c, d, f, BA, 通过 对 推论 1 中 的 
判 据 N < 0 进行 编程 ， 并 利用 Matlab 中 的 线性 矩阵 不 等 式 (LMI) 工具 箱 求解 得 到 ， 在 6 取 不 
同 的 值 时 ， 有 关 时 灌 长 度 +、 判 别 矩 阵 N 的 最 大 特征 值 入 及 Lyaunpov 指数 -ro 的 不 同 数值 ， 
这 些 数 值 被 列 在 表 1 中。 


表 1: 数值 分 析 
0 0 0.01 0.05 0.10 0.15 0.18 0.1835 
T 1.6648 1.5811 1.1287 0.5976 0.3103 0.0519 0.0030 


入 -0.8623 -0.8205 -0.6384 -0.3491 -0.1756 -0.0328 -0.0016 
—ro -0.0908 -1.6733 -0.0941 -0.0847 -0.0674 -0.0234 -0.0013 


从 表 1 中 的 数值 可 以 看 到 ， 随 着 6 的 增 大 ，7 随 之 减 小 ， 入 随 之 增 大 ，-ro 则 先 减 后 增 。 也 
就 是 说 ， 这 个 二 维系 统 的 指数 稳定 性 随 着 灰 参 数 和 矩阵 的 绝对 灰 度 矩阵 的 范 数 5 的 增 大 而 逐渐 减 
弱 ， 且 使 该 系统 保持 指数 稳定 的 绝对 灰 度 矩阵 的 范 数 5 及 相应 的 时 滞 长 度 r 的 最 大 允许 取 值 范 
围 分 别 是 [0, 0.1835] 和 [0.0030, 1.6648]. 


5 ”结论 


通过 对 一 类 中 立 型 灰色 随机 分 布 时 滞 系 统 的 指数 稳定 性 研究 ， 分 别 得 到 了 该 系统 指数 和 鲁 棒 
稳定 的 非 线 性 矩阵 判 据 和 线性 矩阵 判 据 。 针 对 非 线性 矩阵 判 据 ， 提 出 了 一 个 简单 有 效 的 算法 。 
同时 证 明了 该 系统 在 均 方 指数 鲁 棒 稳 定 的 条 件 下 ， 又 是 几乎 必然 指数 鲁 棒 稳 定 的 。 实 例 说 明 ， 
所 给 判 据 是 有 效 的 ， 且 系统 的 指数 稳定 性 和 时 湛 是 随 着 灰 参 数 和 矩阵 的 绝对 灰 度 矩阵 的 范 数 的 增 
加 而 减 小 的 。 该 判 据 在 实际 应 用 方面 是 简单 方便 的 ， 对 区 间 灰 矩阵 分 解 出 的 灰 矩 阵 部 分 ， 不 需 
要 再 将 其 分 解 成 符合 某 种 假设 条 件 的 一 些 抢 阵 的 乘积 形式 ， 只 需要 借助 区 间 灰 矩阵 的 连续 矩阵 
覆盖 集 的 上 下 界 矩 阵 和 一 个 简单 小 程序 就 能 判断 系统 的 指数 稳定 性 。 
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Exponential Robust Stability of Grey Neutral Stochastic Systems 
with Distributed Delays 
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Abstract: In this paper, we investigate the exponential robust stability for a class of grey neutral 
stochastic systems with distributed delays. By constructing the Lyapunov-Krasovskii functional and 
applying the decomposition technique of continuous matrix-covered sets of grey matrix and the It6 
formula, the delay-dependent criteria for exponential robust stability are formulated in the forms of 
non-linear matrix inequalities and linear matrix inequalities, respectively. A generic algorithm for the 
non-linear matrix inequalities is given. The existent criteria for non-linear matrix inequalities are not 
convenient for applications in practice, and our algorithm solves this problem. Numerical examples 
demonstrate the effectiveness of the presented criteria, and that the exponential stability and time- 
delay of systems decrease as the spectral norm of the absolute grey-degree matrix increases. 
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